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1. Introduction 

In the Regge limit, where the total energy s is large while the transfer of 
four-momentum t is low and fixed, the scattering amplitude of two hadrons 
can be rewritten as an exchange of effective particles, i.e. propagating in 
i— channel reggeized gluons, which are also called Reggeons pQ El El El El • 

Due to ordering of generalized leading logarithm approximation contri- 
bution from exchange of A^— Reggeon states is suppressed by factor a^~^ 
where ag is a strong coupling constant. Thus, the leading contribution is 
given by two Reggeon states, i.e. BFKL Pomeron. The equation describ- 
ing this case was firstly derived and solved by Balitsky, Fadin, Kuraev and 
Lipatov [IIIHIE]- Extension of this equation for more than two Reggeons 
was formulated in 1980 [HI EB ^2 by Bartels, Kwiecihski, Praszalowcz and 
Jaroszewicz. It has a structure of Schrodinger equation. The first solutions 
for exchange of three Reggeons appear in nineties jT21 El ^M ^^^ these so- 
lutions correspond to leading contribution to the QCD odderon state as well 
as subleading contribution to the Pomeron state. 

The Schrodinger equation for A^ > 4 Reggeons contains complicated 
colour factor and even in 't Hooft's multi-colour limit j231 ^M ESj, «-e. 
A'^c ^ oo, formulating of its solutions poses real challenge and demands us- 
ing of advanced methods of integrable systems. In the multi-colour limit the 
Reggeon Hamiltonian corresponds to the Hamiltonian SL(2,C) spin magnet 
which turns out to be solvable model. Thus, making use of Q— Baxter and 
Separation of Variable methods one can solve the multi-Reggeon problem 
is multi-colour limit. The solutions for higher A^ = 4, . . . ,8 were obtained 
recently in series of papers ^3 ^Q Ej written in collaboration with S.E. 
Derkachov, G.P. Korchemsky and A.N. Manashov. Similar results for A^ = 4 
appear also in Ref. [T8lll9j . 

The present work is continuation of [201 ^^ ^^11 ^^ C3 ^M El- We 
concentrate on A^ = 4 case with conformal Lorentz spin Uh = 0. Here, 
classification of the four Reggeon states is performed and corrections to the 
WKB calculation are computed. Moreover, we calculate the rich spectrum of 
the the Reggeon energy and the conformal charges {q3,qi} for four reggeized 
gluons. For a broader perspective see also Ref. |21| . 

Thus, in Section 2 we explain when the reggeization of the gluon appears 
|22| . Next, we perform the multi-colour limit (2S10I2S1, discuss the proper- 
ties of the SL(2, C) symmetry and construct invariants of this symmetry, the 
conformal charges. In Section 3 we introduce a Baxter Q— operator method 
[2ni with Baxter equations that allow us to solve the Reggeon system, com- 
pletely. Next, we present the exact solution to the Baxter equations. It con- 
sists in rewriting the Baxter equation into the differential equation, which 
may be solved by a series method. We find the quantization conditions for 
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{53,^4} which come analyticity properties of the Baxter functions. We also 
recapitulate known properties of A^— Reggeon spectrum. The numerical re- 
sults are shown in Section 6 ^HIEI- In particular, we discuss the resemblant 
and winding structures of the 54 and 53 spectrum and also corrections to the 
WKB approximation for ^3 = 0. At the end we make final conclusions. 

2. System with SL(2, C) symmetry 

2.1. Hamiltonian 
In the Regge limit 

s — > 00 and t = const (2.1) 

A^— Reggeon Hamiltonian can be rewritten as a sum of BFKL kernels. Per- 
forming the multi-colour limit ^23, 24, 25], where a number of colours N^ -^ 
00, the system of A^— reggeized gluons is described by Hamiltonian 

Hn = Hn + Hjq, [Hn, Hn\ = (2.2) 

which can be written in terms of the conformal spins ^'1.1\ : 

N N 

Hn = 2_^H{Jk^k+i), Hn = 2_^H{Jk^k+i), (2.3) 

fe=i fc=i 

where 

H{J)=il^{l-J)+i;{J)-2i;{l) (2.4) 

with ip{x) = dlnT{x)/dx being the Euler function and Jn,n+i = Jn,i- Here 
operators, Jfc,fc+i and Jk,k+i, sue defined through the Casimir operators for 
the sum of the spins of the neighbouring Reggeons 

Jk,k+i{Jk,k+i - 1) = (5^') + S(^+i))2 (2.5) 

with Sa = Sa , and Jk,k+i is defined similarly. 

2.2. Symmetry SL{2,C) 

The Hamiltonian (|2.2|1 is invariant under the coordinate transformation 
of the SL(2,C) group 

^, ^ azk + b _, ^ azk+l 

'' czk + d' ^ czk + d 

with k = 1, . . . .,N and ad — he = ad — he = 1. 
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Now, one may associate with the fc— th particle the generators of this 
transformation ^7\ . These generators are a pair of mutually commuting 



holomorphic and anti-holomorphic spin operators, Sa and S^ . They sat- 



isfy the standard commutation relations 






i(-a(3'y5^"' S\ and 



similarly for 5^ . The generators act on the quantum space of the fc-th 
particle, ^(■'fc'^*) as differential operators 



'S'o = Zkd-^^ + Sk , S^_^^ = -d^^ , S'J'' = zp^^ + 2skZk , 

^0 = ZkOzk + Sk , J— = "Cjft J "J-i- = Zj^Ozk "I" ^SkZk , 



where 5j_ = 5} ± iSg while the complex parameters, Sk and Sk, are 
called the complex spins. Thus, the Casimir operator reads 

j^isf^? = (Stf + (4'^si'^ + 5i')5f )/2 = skisk - 1) (2.8) 

and similarly for the anti-holomorphic operator {S )^. 

The eigenstates of SL(2,C) invariant system transform as [281129] 

^{zk,Zk) -^ ^'(^fc,^fc) = {czk + d)-^'^{czk + d)-2^'=^(z^,4). (2.9) 

In statistical physics (|2.Hjl is called the Hamiltonian of the non-compact 
SL(2,C) XXX Heisenberg spin magnets. It describes the nearest neighbour 
interaction between N non-compact SL(2,C) spins attached to the particles 
with periodic boundary conditions. 

For the homogeneous spin chain we have to take Sk = s and Sfc = s. In 
QCD values of (s, s) depend on a chosen scalar product in the space of the 
wave-functions (|2.9|1 and they are usually equal to (0,1) or (0,0) [27|ll9j. 

2.3. Scalar product 

In order to find the high energy behaviour of the scattering amplitude 
we have to solve the Schrodinger equation 

W^=°'^=^)^(zi, Z2, ...,zn) = En'^{zi,Z2, ...,zn) (2.10) 

with the eigenstate ^(zi,Z2, • • • ,zn) being a single-valued function on the 
plane z = (z,^), normalizable with respect to the SL(2,C) invariant scalar 
product 

\m^ = (*|*) = f d^Zid^Z2 . . . d''zN\^{zi,Z2, ...,Zn)W (2.11) 
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where d^Zj = dxidyi = dzidzi/2 with Zi = Zi* . One may notice that it 
is possible to use other scalar products corresponding to different choice of 
{s,s). For farther information see Refs. [dE] and [TH|[Tn]. 

Let us consider the amplitude for the scattering of two colourless objects 
A and B. In the Regge limit, the contribution to the scattering amplitude 
from A^— gluon exchange in the i— channel takes the form 

A{s, t) = is^{ias)^AN{s, t) . (2.12) 

N 

Using the SL(2,C) scalar product (|2.1ip we have 

AN{s,t) = s j d^z^ e-o-P^(|.^(zo)|e-^»^«'v/4||.B(0)), (2.13) 

where the rapidity Y = Ins. Here the Hamiltonian TL;^ is related to the 
sum of A^ BFKL kernels corresponding to nearest neighbour interaction be- 
tween N reggeized gluons. The wave-functions |$^(b)(zo)) = ^A{B){^i — -?o) 
describe the coupling of A'^— gluons to the scattered particles. The zq— 
integration fixes the momentum transfer t = —p . 

2.4- Conformal charges qk and the conformal spins 

The Hamiltonian H2.2|l possesses a complete set of the integrals of motion 
{p,qk} where qk = {Qk,Qk} with k = 2, . . . N are called conformal charges 
while p = {p,p} is the total momentum of the system. In order to construct 
them we introduce the Lax operators jSHl 1^ IHH IHHI in holomorphic and 
anti-holomorphic sectors: 



Lfc(u) = u + i{a-S^''^) 



z5f u - ^si'^ 



'0 ^•5- 

iO^ U — iDg 



Lkiu) = u + ^{a. S^'^) = ( ^ +^^^f _ ^^7 ,,) 1 (2.14) 



with u and u being arbitrary complex parameters called the spectral param- 
eters and (Tq being Pauli matrices. 

To identify the total set of the integrals of motion of the model, one 
constructs the auxiliary holomorphic monodromy matrix 

r^(n)=Li(n)L2(n)...Ljv(n) (2.15) 

and similarly for, the anti-holomorphic monodromy operator T]\[(u). Taking 
the trace of the monodromy matrix we define the auxiliary transfer matrix 
(spectral invariants) 

tNiu) = tr [Tn{u)] = 2u^ + q2U^-'^ + . . . + qN (2.16) 
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and similarly for t]y{u). We see from H2.16p the advantage of using the trans- 
fer matrix tN{u)'- that is a polynomial in u with coefficients given in terms 
of conformal charges qk and g^ , which are expressed as linear combinations 
of the products of k spin operators: 

N 2 

« = -2 E E («i:"*lr' 

«2>«1 = 1 jl=0 

N N 2 

- _ V^ \^ \^ c- f oih) Q(i2)\ { oiis) qM 

94 - Z^ Z^ 2-^ £ili2i3i4 [^j, •3ji J \^^J2 ^J2 

*2>«1=1 «4>J3 = 1 il,i2=0 

N N N 2 

96 = -^2^ Z^ 2^ Z^ ^JlJ2J3i4J5i6 (^'5^^i -S^ji j 

■i2>«l = l«4>«3 = l«6>«5 = lil,J2,i3=0 

where £4^12. ..i^ is completely anti-symmetric tensor and eiiJ2--«fc ~ 1 fo'^ ^1 < 
12 < ■ ■ ■ < ik- So for even k we have a formula for conformal charges 

N 2 k/2 

^2 > «1 = 1 JiJ2'-.ifc/2=0 n=l 

U > «3 = 1 

^n ^ ^n— 1 — J- 



:<(n) ofe) o(«3) 



(2.18) 
For odd fc's we have 

;5 _ 9 V^ c- f qin) qin) qi-i-s) 

ya — ^ /__^ tn«2«3 I '-'0 '-'1 '-'2 

*lj«2,i3 = l 

TV 2 

- 2 V \^ r- ■ ■ F- ■ ■ f Q^^i) 9^*2) Q^*' 

~ ^ Z_^ Z_^ tn«2«3tjlJ2J3 \^^ji '-'j2 '^is 

is >i2 >n =1 ii j2 ,i3=o 

Af N 2 

- _ 9 \^ \^ \^ =- /'c<(n)c'fe)c.(i3)\ /'c.(i4)c(J5) 

«3,J2,n=l i5>«4=l i4=0 

TV AT AT 2 

A, - "S^ \^ V V P ( qin) q(i2) qii3)\ 

97 — ^^ /_^ 2_^ 2_^ tj^j2«3«4«5«6«7 \^'-'0 '-'1 '-'2 j 

«3,«2,«1 = 1 J5>«4 = 1 i7>«6 = l i4 J5=0 

X (•5j:*^s];^^) (5]:«^4:^)) (2.19) 
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and the general expression for an odd number of the conformal spins is 
^ _ 2(-l)(^+i)/^ ^ ^ 2 

V 2 /• i3,«2,n=l Z5 > Z4 = 1 iiJ2v,i{fc-3)/2=0 

Z7 > Z6 = 1 

ik > 4-1 = 1 

(fc-3)/2 
^(n)^(^2)^(^3)^ -Q J^^fe„+2)^ji2„+3)^ _ (2.20) 

n=l 

In the above formulae we have two basic blocks (Sq 5}*^ 5*2 ) and 
■^ii 71 ) whose products are summed with antisymmetric tensor eiiJ2...ij.. 

2.5. Two-dimensional Lorentz spin and the scaling dimension 

The Hamiltonian H2.3p is a function of two-particle Casimir operators 
|27| . and therefore, it commutes with the operators of the total spin Sa = 

^^ Sa and Sa = Ylk '^a ) acting on the quantum space of the system 
Vn = y(^i'^i) y(«2,s2) ... y(siv,^iv). This imphes that the eigenstates 
can be classified according to the irreducible representations of the SL(2, C) 
group, yC^'^), parameterized by spins {h,h) j2'7j . 

The Hamiltonian depends on differences of particle coordinates so the 
eigenfunctions can be written as 

^p-(zi, Z2, ...,zn)= / d'^zoe"°-P'^{zi - zo, Z2 - zo, . . . ,zn - zq) . (2.21) 

The eigenstates '^{zi,Z2, . . . ,zn) belonging to V^^'^' are labelled by the 
centre-of-mass coordinate zq and can be chosen to have the following the 
SL(2, C) transformation properties 

^{{z,'-z,'}) = 

= (czo + df'iczo + df'^ (fliczk + df'^icz, + df^A ^({4 - ^-o}) 

(2.22) 

with Zo and 'zq transforming in the same way as Zk and Zfc, H2.6p . As a 
consequence, they diagonalize the Casimir operators: 

{S^ - h{h - l))*(zi, Z2, . . . , zjv) = (2.23) 
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corresponding to the total spin of the system, 

N 2 N 

s"=Y.Y. t^sf'^ = -^2 - E ^^(^^ - 1) • (2-24) 

«2,*i=ii=o fc=i 

The complex parameters {sk,'Sk) and {h, h) parameterize the irreducible 
representations of the SL(2, C) group. For the principal series representation 
they satisfy the conditions 

Sk-Sk = ris^ , Sk + (sk)* = 1 (2.25) 

and have the following form 

1 + Us, _ 1 — Ug, 

Sk = Y^ + ivs^ , Sk = Y^ + ivs^ (2.26) 

with 1/sj.being real and n^j. being integer or half-integer. The spins (/i, h) 
are given by similar expressions with nsj.and Ugk replaced by Uh and Vh, 
respectively 

h = — ^ ViVh, h = — ^ ViVh- (2.27) 

The parameter Ug^ has the meaning of the two-dimensional Lorentz spin of 
the particle, whereas Ugk defines its scaling dimension. To see this one can 
perform a 27r-rotation of the particle on the plane, and find from eigenstates 
transformations H2.22|l that the wave-function acquires a phase. Indeed 

^fc ^ ^fc e^^* and Zfc ^ Zfc e-^'^* gives ^(z^, z^) ^ (-l)2"=fc*(zfc, z^). 

(2.28) 
For half-integer n^j. it changes the sign and the corresponding representation 
is spinorial. Similarly, to define scaling dimension, s + 1 = 1 + 2ivsf, one 
performs the transformation 

z ^ \z and 'z ^ Xz giving ^(zfc,^fe) ^ A^"''^*'''''=^(2;fc, z^) . 

(2.29) 
Because the scalar product for the wave- functions is invariant under SL(2, C) 
transformations, (|2.6p . the parameter Vsk is real. 

We notice that the holomorphic and anti-holomorphic spin generators as 
well as Casimir operators H2.5p are conjugated to each other with respect to 
the scalar product H2.1ip : 

[5(^)]t = _5(f), [j,]t = i_j^. (2.30) 
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Moreover, because of the transformation law (j2.:-{()|l . h* = 1 — h ^. This 
imphes that W^ = H^ and, as a consequence, the Hamiltonian is hermitian 
on the space of the functions endowed with the SL(2,C) scalar product, 

TCj^ = TCn- 



2.6. Conformal charges qk as a differential operators 

We noticed in the previous Section that the conformal charge operators 
Qk are given by invariant sum of linear combinations of the products of k 
spin operators. They can be rewritten as k-th order differential operators 
acting on (anti)holomorphic coordinates {z,z). 

Two particle spin square can be written as 

31=0 

= -2^^n - Zi2fdz,^d;,^^ + (zii - Zi^){si^d^^^ + Si^d^^^) + siS2 . (2.31) 

For homogeneous spins s = si = S2 = ■ ■ ■ = sn, what is also the QCD case, 
we have 

TV / 2 

*2>n=i \ji=o 
N 
= E {iz^,^lf^'-'^^.^^^.^^{z.,,,r + 2s{s-l)) , 
i2>«i=l 
N 
/; _ 9 V^ c- C'(«i)c<(«2) c-(J3) 

«1,«2,«3 = 1 

N 

i3>i2>n=l 
+ SZJ2J3(2J3Ji — ^1-^12) Oz^Ozi + SZi^j^i^yZi^j^i^ — Zi^i^jOzi Ozi 



,2^ o o„2^ o o„2. 



■^■^ Zi-^i2C'zi ^S Zi2i:jC'z^ ^S Zi^i-^Oz^ I , (Z.6Z) 

where Zij = Zi — Zj. Similar relations hold for the anti-holomorphic sector. 



- denotes complex conjugation 
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In that way one can also construct operators for the higher conformal 
charges. They have a particularly simple form for the SL(2,C) spins s = 

% = ^ /_^ ^jij2 ■ ■ ■ ^jk-i,jk^jk,ji'^Z]i ■ ■ ■ '^Zj^_^c)zj^ (2.33) 

i<ii<J2<---<ife<Af 

as well as for s = 1 

qk=i'' Yl ^'n ■ ■ ■ ^^Jk-i ^^3k ^hh ■ ■ ■ Zj,_i,h ^h ,h ■ (2-34) 

i<h<J2<...<jk<N 

The eigenvalues of the lowest conformal charge q2 can be parameterized by 
a conformal weight h H2.27|l and complex spin s (|2.26p as follows 

q^ = -h{h-l)+Ns{s-l) (2.35) 

2. 7. Other symmetries 
The states H2.2H1 have additional symmetries |27| : 

'^^q,q{zi,Z2,...,ZN) = ^g,g(z2, Z3, . . . , Zi) 

= eie^(M)^^_(^^,^-,...,^lv), 

M*±(zi,Z2,...,^iv) = *^(Z7V,^7V-I,...,i'l) 

= ±^±(zi,Z2,...,i7v) (2.36) 

so called cyclic and mirror permutation where the conformal charges are 
denoted by g = (^'2,93, • • • ,9n) and q = (^2) 93) • • • j^n)- The generators 
P and M, respectively, commute with the Hamiltonian TL but they do not 
commute with each other. They satisfy relations 

P^ = M^ = 1, pt=p-i=p^-i, M'^=M, PM = MP"^ =MP^~^ 

(2.37) 
The phase 6n{q) which is connected with eigenvalues of P is called quasi- 
momentum. It takes the following values 

eN{q,q) =27r^, foi k = 0,1, . . . ,N - 1 . (2.38) 

The eigenstates of the conformal charges qk diagonalize Ti. and P. 

The transfer matrices (|2.16|1 are invariant under the cyclic permutations 
PT tN^u) P = tiy{u) whereas they transform under the mirror transformation 
as 

Mtjv(n) M = (-l)^t7v(-n). (2.39) 
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Substituting (|2.1f)jl into H2.,'-if)|l one derives a transformation law of the con- 
formal charges qk 

¥Uk^ = qk MqkM={-lfqk (2.40) 

and similarly for the anti-holomorphic charges. Since the Hamiltonian (|2.,S|I 
is invariant under the mirror permutation, it has to satisfy 

W(gfcifc) = MW(gfeifc)M = W(MgfcM,MtfcM) = n{{-lfqk, {-ifW ■ 

(2.41) 
This implies that the eigenstates of the Hamiltonian H2.Hp corresponding to 
two different sets of the quantum number {qk^Qk) ^^^^ {( — 1)'^^^, ( — 1)^^;,} 
have the same energy 

EN{qk,qk) = EN{{-lfqk, {-ifqk) ■ (2-42) 

Similarly one can derive relation for quasimomentum 

ON{qk,qk) = -ON{{-l)''qk, (-1)%) • (2.43) 

The cyclic and mirror permutation symmetries come from the Bose sym- 
metry and they appear after performing the multi-colour limit |2,S| . Physical 
states should possess both symmetries. 

3. Baxter Q-operator 

3.1. Definition of the Baxter Q-operator 

The Schrodinger equation (|2.1()|l may be solved applying the powerful 
method of the Baxter Q-operator ^^. It depends on two complex spec- 
tral parameters u, u and in the following will be denoted as Q{u,u). This 
operator has to satisfy the following relations 

• Commutativity: 

[Q(n,n),Q(t;,U)]=0, (3.1) 

• Q — t relation: 

[iNiu),q{u,u)] = [f^(u),Q(n,n)] =0, (3.2) 

• Baxter equation: 

iNiu)Q{u,u) = {u + is)^Q{u + i,u) + {u- is)'^Q{u - i,u) , (3.3) 

iNiu)Q{u,u) = {u + is)'^(Q{u,u + i) + {u - i-s)^Q{u,u - i) , (3.4) 
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where tAf(ti) and t]y{u) are the auxihary transfer matrices (j2.16p . According 
to H3.2|l the Baxter Q(ti,Il)-operator and the auxihary transfer matrices as 
well as the Hamiltonian H2.2p share the common set of the eigenfunctions 

Q{u,u)'i>g^g{zi, Z2, . . . , Zn) = Qq;q{u,u)^q;q{zi,Z2,...,ZN). (3.5) 

The eigenvalues of the Q-operator satisfy the same Baxter equation H,S.,S|1 and 
H3.4p with the auxiliary transfer matrices replaced by their corresponding 
eigenvalues. 

In the paper J2jj the Q-operator was constructed as an A^— fold integral 
operator 

Q(n,M)^(zi,Z2,--- iZm) = I (fwi / <f'W2 ... (fwN 

Qu,u{^l,Z2,... ,Zn\wi,W2,... ,WN)^iwi,W2,...,WN) , (3.6) 

where the integrations are performed over two-dimensional Wi— planes. The 
integral kernel in (|3.6|1 takes two different forms: 



QS(-» 



N r il-2s 

a{2 — 2s,s + iu,s — lu) n I I 



^^ , nl — s— m r il — s+m 

k=i [Wk - Zk\ [Wk - Zfc+lJ 

(3.7) 
and 

oa(4»T = n -, — '^t^^'^"' ,-.. . (3^8) 

k=i [Zk - Wk\ [Zk - Wk+l\ 

which appear to be equivalent 22I- In Eqs. H3.7p and (|3.8p function a(. . .) 

factorizes as 

-p /-| — \ 

a{a, j3, . . .) = a{a)a{(3) . . . and a{a) = — — (3.9) 

r(a) 

and a is an anti-holomorphic partner of a satisfying a — a G Z. Moreover, 
the two-dimensional propagators are defined as 

[zk - Wk]'"" = {zk - Wk)'"' {zk - Wk)'"' ■ (3.10) 

In order for the Baxter Q— operators to be well defined, (|3.7p and H3.8p . 
should be single- valued functions. In this way we can find that the spectral 
parameters u and u have to satisfy the condition 

i{u — u)=n (3-11) 
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with n being an integer. 

The Baxter Q-operator has a well defined pole structure. For Q^+^ (u, u) 
we have an infinite set of poles of the order not higher than A^ situated at 



{ 



^)}; { 



u„ 



—i{s — m),u^ = —i{s — m)} 

(3.12) 

with 7n,m = 1,2,.... The behaviour of Qqq{u,u) = Q^-'{u,u) i.e. an 

eigenvalue of Q^'^'{u,u), in the vicinity of the pole at m = fn = 1 can be 
parameterized as 



^± 



±1 „ ^\ 



Qq,q{ui + e, Ui +e) = R {q, q) 



7N + 



-E^{q,q) 



+ ... 



(3.13) 



The functions R^{q,q) fix an overall normalization of the Baxter operator, 
while the residue functions E^{q,q) define the energy of the system (see 
Eqs. H3.21|) and H3.24|l below). It has also specified asymptotic behaviour. 
For |ImA| < 1/2 and ReA -^ cx) 

Qgq{X-in/2, X+in/2) ~ e^®'^^'?'^) x^+h-Nis-s) ^^-ieUg,Q) ^i-h+i-h-N(s-s) ^ 

(3.14) 
where @h is a phase that should not be confused with quasimomentum 

3.2. Observables 
The Hamiltonian H2.2p may be written in terms of the Baxter Q-operator 



"Hn = eAT + Z-— In Q*-"*"-* {u + is,u + is] 
du 



t-— in 
du 



M=0 



U — tS,U — IS) 



u=0 



where the additive normalization constant is given as 

ejv = 2NRe['tp{2s) + V(2 - 2s) - 2V'(1)] . 
Applying to l|3.15p the eigenstate "ifg we obtain the energy 



d 



EN{q, q) = sn + i-f- In [Qg^giu + is,u + is) {Qq^q{ 



du 



U — IS, U — IS 



)y 



(3.15) 



(3.16) 



u=0 

(3.17) 
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or equivalently 

J r 

u'^'^Qq,g{u + i(l -S),U + i{l - S)) 
X Q_q^_q{u + i(l - S),U + i(l - s)) 



EN{q,q) = -Im-— In 
du 



(3.18) 



M=0 



where Qqq{u,u) = Qgg{u,u) is eigenvalue of the Q^~^'{u,u) operator, while 



!q,q\^, "■; — ^qq 

±q = {q2,±q3,...,{±fqN) (3.19) 

are the conformal charges. 

It is also possible to rewrite the quasimomentum operator in terms of 

Moreover, using the mirror permutation (|2.40p one finds the following parity 
relations for the residue functions R~^{q,q) defined in H3.13|) : 

R+{q,q)/R+{-q, -q) = e^^^^^^'^) (3.21) 

and for the eigenvalues of the Baxter operator: 

Qq,q{-U, -U) = e^^^('?'^) Q-q,-q{u, u) , (3.22) 

where —q = (52, — ^3, . . . , ( — l)"g'n) and similarly for q. Examining the be- 
haviour of H3.22|l around the pole at u = n^ and u = Ui and making use of 
Eq. (|3?T3|l one gets 

R^{q,q) = (-l)^e*«-(«'5)^T(_^^_^)^ ^±(g^^) = _i^T(_^^_^). 

(3.23) 
To obtain the expression for the energy Efq{q,'q), we apply H3.17p and 
replace the function Qq^q{u ± i{l — s),u± i{l — s)) by its pole expansion 
H3.13J) . Then, applying the second relation in (|3.23|l . one finds 

EN{q,q) = E+{-q, -q) + {E+ {q,q)Y = Re [E+{-q, -q) + E+{q,q)] , 

(3.24) 
where the last relation follows from hermiticity of the Hamiltonian H2.3|l . 
We conclude from Eqs. I|3.24|l and (|3.13|l . that in order to find the energy 
E]\f{q,q), one has to calculate the residue of Qq^q(u,u) at the {N — l)-th 
order pole at u = i(s — 1) and u = i{s — 1). 
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3.3. Construction of the eigenfunction 

The Hamiltonian eigenstate ^p,{q,q}iz) is a common eigenstate of the 
total set of the integrals of motion, p and {q,q} as well as the Baxter Q- 
operator. Thanks to the method of the Separation of Variables (SoV) de- 
veloped by Sklyanin (HHl l2lj we can write the eigenstate using separated 
coordinates x = {xi, . . . , afjv-i) as 



x($q,g(fi,...,f^_i))* , (3.25) 
where U^^ is the kernel of the unitary operator while 

{^q,q{xi,...,XN-l))* = 



N-l / -n/ , ■ \-n/- — ^ \ N 



giejv{g,5)/2 



^^^[nr^7TJ^^;;:jnr^j^W)) ^-^(---^)- (3-26) 

The functions Qq,q{xk,Xk) are eigenstates of the Baxter Q-operator. In 
contrast to the Zi = {zi,^i) - coordinates, the allowed values of separated 
coordinates are 

ink _ iiT-k 

Xk = Vk- — , Xk = i^k + — (3.27) 

with Uk integer and Uk real. Integration over the space of separated variables 
implies summation over integer n^ and integration over continuous I'k 






N-l / oo ,_, 

^ _ ^ |2 
Xk Xj I , 

fc=l \nf,=-oo ■" -°° / ^" '^' j,k=l 

j>k 

(3.28) 
where \xk - Xj\ = {uk - Uj)"^ + {rih - nj)"^ /A. 
The integral kernel f/-.^ can be written as 

[/.^(zi, ...,zn)= cn{x){p^)^''-^^'^ j cfwNC^'P-'^^U^izi, ...,zn; wn) , 

(3.29) 
where 2p- wn = PWN + pwN, 

U^{zi, ...,Zn; wn) = ^M-\.(^.AN-2.(g.,) ■ ■ ■ K.(L_.) (^1' • • • ' ^Mwn) 

(3.30) 



'^N-l,ixiy^N-2,{x2) ■ ■ ■ ''■l,(a;jv_i) 
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for even N, and 

TT t^ - -> \ TaC^is) . (1-s,1 — s) . (1-s,1-s)1 f^ _ I ^ N 

U^{zu..., zn;wn) = [^N-i,(x^)^N-2,ix2) ■ ■ ■ KixM-i) \ (^1' ■ ■ ■ ' zn\wn) 

(3.31) 
for odd A'^. Here the convolution involves the product of {N — 1) functions 
^N-k,{xk) with alternating spins (s, s) and (1 — s, 1 — s). They are defined 
as 



A' 



(s,s) j -t -> \-' -• \ r 1— x+i« 

7V-n,(x)(^"' • • • ' ^A^l^n+l' ■ ■ ■ .Vn) = Vl - ^2] 

X I JJ \zk- Vk] "" ™ [^fe - Vk+i] ^^™ I [^AT - Vn] ^ ™ , (3.32) 

\k=n+l / 

where the convolution [Ajv-fc,(xfc)-'^7V-fc+i,(xfc_i)] contains (N — k) two-dimen- 
sional integrals. The coefficient cn{x) is given for A^ > 3 

[iN-l)/2] 

cn{x) = JJ ia{s + ix2k,s - ix2k)) 
k=l 

[N/2-1] 

X JJ (a(s-Fix2fc+i,s-ix2fc+i)) (3.33) 
fc=i 

while the products go over integer numbers lower than upper limit. For 
A^ = 2 we have 02(^1) = 1. 



4. Quantization conditions in the Q-Baxter method 

In Ref. 23 the authors describe a construction of the solution to the 
Baxter equations, (|3.3p and (|3.4p . which satisfies additionally the conditions 
H3.12J1 and H3.14|l . This can be done by means of the following integral 
representation for Qq^q{u,u) 

Qq,-q{u,u)= /'^z-»z--Q(z,z), (4.1) 

J ^^ 

where we integrate over the two-dimensional z— plane with 'z = z* and 
Q{z,z) depends on {q,q}. The advantages of this ansatz are: 

• the functional Baxter equation on Qq^q{u,u) is transformed into the 
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A'^— th order differential equation for the function Q{z,z) 
'z'izd,fz^-' + z-'{zd,fz'-^ 



N 



-2izd,f -^i^Qkizd, 



k=2 



Q{z,z)=0. (4.2) 



A similar equation holds in the anti-holomorphic sector with s and qk 
replaced by s = 1 — s* and Qk = ll^ respectively. 

• the condition H3.11|l is automatically satisfied since the z— integral in 
the r.h.s. of (|4.ip is well-defined only for i{u — u)=n. 

• the remaining two conditions for the analytical properties and asymp- 
totic behaviour ofQq^q{u,u), Eqs. (|,S.12p and (|,S.14|1 . become equivalent 
to the requirement for Q{z,z = z*) to be a single- valued function on 
the complex z— plane. 

The differential equation H4.2p is of Puchsian type. It possesses three 
regular singular points located at z = 0, z = 1 and z = oo. Moreover, it has 
N linearly independent solutions, Qa{z). The anti-holomorphic equation 
has also A^ independent solutions, Qbi'z). 

Now, we construct the general expression for the function Q{z,'z) as 

N 

Qiz,z)=Y,Qaiz)CabQi,iz), (4.3) 

a, 6=1 

where Cab is an arbitrary mixing matrix. The functions Qa{z) and Qb{z) 
have a nontrivial monodromy^ around three singular points, z, z = 0, 1 and 
oo. In order to be well-defined on the whole plane, functions Q{z,'z = z*) 
should be single-valued and their monodromy should cancel in the r.h.s. 
of l|4.3p . This condition allows us to determine the values of the mixing 
coefficients. Cab, and also to calculate the quantized values of the conformal 
charges qk- 

The differential equation H4.2|l is also symmetric under the transforma- 
tion z -^ l/z and q^ — > (— l)'^^^. This property is related to Eq. (|3.22p and 
leads to 

Qq,q{z,z) = e*^^('''5) Q_^^_^(l/z, l/z) , (4.4) 



^ The monodromy matrix around 2: = is defined as Qn (ze^^') ~ M„kQl. (2) and 
similarly for the other singular points. 
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where ±q = (^2, ±93, • • • , (±)''^'?Ar) denotes the integrals of motion corre- 
sponding to the function Q{z,z). The above formula allows us to define the 
solution Q{z,z) around z = co from the solution at z = 0. Thus, applying 
H4.4p we are able to find Q{z,z) and analytically continue it to the whole 
z— plane. 

4-1- Solution around z = 

We find a solution Q{z) ~ z" by the series method. The indicial equation 
for the solution of Eq. H4.2J1 around z = reads as follows 

(a-l + s)^ = (4.5) 

and the solution, o = 1 — s is A^— fold degenerate. This leads to terms 
~ Log (z) with k < N — 1. We define the fundamental set of linearly 
independent solutions to (|4.2p around z = as 



Q^l'^z) = z'-'miz), 

Qt\z) = z'-^[m{z)Log"^~\z) 



m~l 



+ Y.cl_,uu+i{z)hog^-^-\z) 



k=l 

with 2 < m < N and where for the later convenience 



(4.6) 



k _ {m - 1)! 
"""^-^ ~ (klim - k - ly.)' ^ ' 

The functions Um{z) are defined inside the region \z\ < 1 and have a form 

«„,(z) = l + 5;z"n(r)(g). (4.8) 

n=l 

Inserting H4.6p and H4.8|l into (|4.2|1 . one derives recurrence relations for 

Un (q)- However, in order to save space, we do not show here their ex- 
plicit form. 

In the anti-holomorphic sector the fundamental set of solutions can be 
obtained from (|4.6|1 by substituting s and g^ by s = 1 — s* and q^ = qI, re- 
spectively. Sewing the two sectors we obtain the general solution for Q{z,'z) 
around z = as 

N 

Q{z,z)^'^=' Y. Q^^H^)c!SnQm{z). (4.9) 

m,m=l 
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The above solution H4.9|l should be single- valued on the z— plane. Thus, im- 
posing single- valuedness condition on (|4.9p we find a structure of the mixing 



matrix C^^ which for n + ?n, < A^ + 1 



C(o) 



a 



N —n—m+l 



(n-l)!(m-l)! 



fc=0 



k\ 



■ Oik+n+rn—l 



(4.10) 



with 0", Oi, . . . ,aAr_i being arbitrary complex parameters and a^r = 1. Be- 
low the main anti-diagonal, that is for n -|- TTi > A^ -|- 1, Cnm vanish. 

The mixing matrix C^:^ depends on A^ arbitrary complex parameters 
a and Uk- However, two parity relations, Eqs. (|H.2ip and (|4.4p . fix ct = 
ex.p{i6iy{q,q)), with 9n{q,q) being the quasimomentum, and lead to the 
quantization of the quasimomentum. Later, we will use l|4.4p to calculate 
the eigenvalues of 0]^{q,q) (see Eq. (|4.25p 'l. 

The leading asymptotic behaviour of Q{z,'z) for z — > can be obtained 
by substituting (|4.inp and H4.6p into (|¥^ . It has a form 



Qq,q{z,z)=Z "Z 



1-s-l-s ^iOr^{q,q) 



Log^-^(2z) Log 

(A^-1)! 



N-2, 



ZZ] 



■aN-i+ 



hogizz) 



1! 



(A^-2)! 

(l + 0(z,z)). (4.11) 



Making use of the integral identity 



d^z • 

\z\<p ^^ 



tt5. 



m—s—iu.m- 



-l)"n! 



(m 



lU 



\n+l 



+ 0{{m- s 



lU) 



(4.12) 



with m and m positive integer, we can calculate the contribution of the 
small— z region to the eigenvalue of the Baxter equation H4.ip . The function 
Qq^q{u, u) has poles of the order A^ in the points u = i{s—m) and u = i(s—fn) 

1 one finds from (14.1111 



what agrees with H3.12p . For m 

Qq,q{ut + e,ut + e) = 



m 



(ie) 



N 



[1 + ie aN-i + ... + {ief'^ Q2 + {ief-' ai + 0(e^)] 



(4.13) 



where uf and u^ are defined in (|3.12p . One can see that the integration in 
H4.ip over the region of large z with H4.4p and (|4.1ip gives the second set of 
poles for Qq^g{u,u) located at u = — i(s — m) and u = —i{s — m). 
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Comparing H4.1H|1 with (|H.i:-{|l one obtains 

R-^{q,q) = -^e^^-(^'5) , E+{q,q) = aN-iiq,q) • (4.14) 

Now, we may derive expression for the energy 

EN{q,q) = Re [aN-i{-q, -q) + aN^i{q,q)] . (4.15) 

The arbitrary complex parameters On, defined in H4.1())l . will be fixed by the 
quantization conditions below. 

In this Section we have obtained following Ref. ^\ the expression for 
the energy spectrum E]\[{q,q), as a function of the matrix elements of the 
mixing matrix H4.1()|l in the fundamental basis (|4.6|1 . Moreover, we have 
defined the solution to the Baxter equation Q{u,u) and reproduced the 



analytical properties of the eigenvalues of the Baxter operator, Eq. H3.12p . 

4-2. Solution around z = 1 

Looking for a solution of (|4.2|1 around z = 1 in a form Q{z) ~ (z — 1)^ 
we obtain the following indicial equation 

7V-3 

(b + 1 + h - Ns){b + 2 - h - Ns) Ylib - k) = , (4.16) 

fc=0 

where h is the total SL(2,C) spin defined in H2.35J1 . Although the solutions 
b = k with k = 0, . . . ,N — 3 differ from each other by an integer, for 
/i 7^ (1 + n/j)/2, no logarithmic terms appear. The Log(2:)— terms are only 
needed for Im/i = where the additional degeneration occurs. 

Thus, we define the fundamental set of solutions to Eq. I|4.2jl around 
z = 1. For Im/i / it has the form 

q51)(z) = z^-'{l - zf'^^'~\i{z) , 

Q^^\z)=z'-^{l-zf^^^-'v,{z), 

QW{z)=z'-'{l-zr-\^{z), (4.17) 

with m = 3, . . . ,N. The functions Vi{z) [i = 1,2) and Vm{z) given by the 
power series 

CO oo 

v.{z) = i + Y,{^-zrv^\q), vm{z) = i+ Y. (i-^r4™Hg), 

n=l n=N—m+l 

(4.18) 
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which converge inside the region |1 — z| < 1 and where the expansion coef- 
ficients Vn and Vn satisfy the A^— term recurrence relations^ with respect 
to the index n. For /i = (1 + nh)/2 G 2Z + 1, one Log(z)— terms appear so 
for Uh > 0, 

h={l+nh)/2 

= zi-*(l-z)^'^-("'^+3)/2[(l-z)"hLog(l-z)t>2(z)+iyi(z)] , (4.19) 

where the function V2{z) is the same as before, vi{z) = Yl'k=o^kZ^ and 
the coefficients v^ satisfy recurrence relations with the boundary condition 
Vnt^ = 1- For /i G Z we have two additional terms: Log(2;) and Log^(2:). 

Similar calculations have to be performed in the anti-holomorphic sector 
with s and h replaced by s = 1 — s* and h = 1 — h* , respectively. A general 
solution for Q{z,'z) for Im(/i) / with respect to the single- valuedness can 
be constructed as 

Q(z, z) ^'t' (ihQ^^\z)Q^l\^) + P,.hQ^^\z)Q^^\z) 

N 

+ E Q^rkHz)lmmQS{-^)- (4.20) 

m,rn=3 

Here the parameters Ph and 7mm build the C^^-* matrix where Q{z,z) = 
Qm CmmQm ■ The /3— coefficients depend, in general, on the total spin h 
(and h = 1 — h*). They are chosen in H4.2()|l in such a way that the symmetry 
of the eigenvalues of the Baxter operator under h ^ 1 — h becomes manifest. 
Thus, the mixing matrix C^^' defined in (|4.2()|1 depends on 2-|-(A^ — 2)^ com- 
plex parameters /3/j, Pi-h and 7mm which are some functions of the integrals 
of motion {q,q), so, they can be fixed by the quantization conditions. 

For h = {l+nh)/2 the first two terms in the r.h.s. of H4.2()|l look differently 



in virtue of M.l 

Q{z,z 



= /3i Qriz)Qriz) + Q','\z)Q\'>{z) 

h=il+ni,)/2 



+ /32Q^'^(z)Q^'\z) + ..., (4.21) 



where ellipses denote the remaining terms^. 



^ The factor z^~" was included in the r.h.s. of 14.1711 and Il4.1t)|l to simplify the form of 

the recurrence relations. 
* Equation 14.2111 describes solutions only for h = [1 + nh)/2 where Uh G 2Z. 
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Substituting (|4.2()|l and (|4.21|1 into H4.HI and performing integration over 
the region of |1 — z| <C 1, one can find the asymptotic behaviour of Q{u,u) 
at large u. 

Let us consider the duahty relation H4.4p . Using the function Q{z,^) we 
evaluate H4.2()p in the limit \z\ — > 1. In this way, we obtain set of relations for 
the functions Pi{q,q) and 'ymm{Q,Q)- The derivation is based on the following 
property 

N 

Q^^\l/z--q)=Y,SahQ^^\z-^q). (4.22) 

6=1 

with Im(l/z) > and where the dependence on the integrals of motion 
was explicitly indicated. Here taking limit z — > 1 in (j4.17p and H4.18|l and 
substituting them to (|4.22p we are able to evaluate the 5— matrix 



Q „—iTT(Ns — h—l) n „—iw(Ns+h—2) n / i \fc— 3 

■^11 = e , J22 = e ' , bk^k+m = i-J^j 



(A;-2s-l), 



ml 

(4.23) 

with {x)m = r(x + m)/r(x), 3</i;<iVandO<m<Af-A;. Similar 

relations hold in the anti-holomorphic sector, 

"C JiriN-s-h-l) -Q jTT(Ns+h-2) "c / 1^fc-3 (^ — 2s — l)m 

ml 

(4.24) 

The S*— matrix does not depend on z because the Q— functions on the both 

sides of relation H4.22p satisfy the same differential equation l|4.2p . 

Now, substituting l|T^ and (|0^ into !f^ . we find 

Ph{q,q) = e^''-('?'5)(-l)^-»+-'^/3,(-(7,-g), 

Af 
ln,rn{q,q) = e*^^^'?'^) Yl Snmlnn{-Q,-q)Srm- (4.25) 

n,n>3 

In this way, similarly to the energy, Eq. (|4.15p . which was calculated from 
the mixing matrix at z = 0, the eigenvalues of the quasimomentum, 0]\f{q,q), 
can be calculated from the mixing matrix at z = 1, from the first relation in 
H4.25p . In the special case when q2k+i = 92fc+i =0(^ = 1, 2...), this means 
Ph{q,q) = Ph{—q, —q), the quasimomentum is equal to 

Qi0N{q,q) — l^_-^\Nns+nh _ u 2Q) 

4-3. Transition matrices 

In the previous Sections we constructed the solutions Q{z,^) to (|4.2p 
in the vicinity of z = and z = 1. Now, we glue these solutions inside 
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the region |1 — z] < 1, \z\ < 1 and, then analytically continue the resulting 
expression for Q{z,^) into the whole complex z— plane by making use of the 
duality relation H4.4p . 

Firstly, we define the transition matrices il(g) and ^{q): 

N N 

Q^f) (z) = ^ n„,^{q) Q« {z) , Qi°) (z) = Y, ^nm{q) Q^^^ {z) . (4.27) 

which are uniquely fixed ^3- The resulting expressions for the matrices 0(g) 
and r2(g) take the form of infinite series in q and g, respectively. Substituting 
H4.27|l into (|4.9p and matching the result into (|4.2n|l . we find the following 
relation 

&^\q,q) = [a{q)fC^^\q,q)Tl{q). (4.28) 

The above matrix equation allows us to determine the matrices C^^> and 
C^^> and provides the quantization conditions for the integrals of motion, q^ 
and g^j with k = 2>, . . . ,N . Therefore, we can evaluate the eigenvalues of the 
Baxter Q— operator, Eq. (|4.H1 . Formula (|4.28p contains A^^ equations with: 

• {N — 1) a— parameters inside the matrix C^^\ 

• 2 + (A^ — 2)^ parameters (5i^2 and 7m,m inside the matrix C^^' , 

• {N — 2) integrals of motion q^, . . . ,qN where q^. = ql 

Thus, we obtain (2A^ — 3) nontrivial consistency conditions. 

The solutions to the quantization conditions (|4.28p will be presented in 
detail in next Sections. 

5. Properties of the AT— Reggeon states 

In this Section we characterize the spectra of the conformal charges 
obtained by numerical calculations ^3IH1|- Here, parametrization of the 
spectra is presented and the spectrum symmetries are shown. Moreover, 
descendent states are described. 

5.1. Trajectories 

Solving quantization conditions H4.28p we obtain continuous trajectories 
in the space of conformal charges. They are built of points, {q2{i'h), ■ ■ ■ , QNi'^h)) 
which satisfy H4.28p and depend on a continuous real parameter v^ entering 
(Z2, (|2..S5p and H2.27p . In order to label the trajectories we introduce the set 
of the integers 

i = {iui2,...,l2iN-2)} (5.1) 
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which parameterize one specified point on each trajectory for given h. Spe- 
cific examples in the following sections will further clarify this point. 

Next we calculate the observables along these trajectories, namely the 
energy (|4.15p and the quasimomentum H4.25|l . The quasimomentum is con- 
stant (|2.38p for all points situated on a given trajectory. The minimum of 
the energy, which means the maximal intercept, for almost all trajectories is 
located at Uh = 0. It turns out that the energy behaves around I'h = like 

ENii^h-, ^^'•°"°^) = E%°'''''' + aNUh^ + 0{uh^) (5.2) 

Thus, the ground state along its trajectory is gapless and the leading contri- 
bution to the scattering amplitude around Uh may be rewritten as a series 
in the strong coupling constant: 

A{S, t) ~ -is Y, (^«s)^7= 7-^JJ^U,NmB,Nit) , (5.3) 



7V=2 



(asO-Arlns)V2 



where Og = OsNc/tt and ^x,N{t) are the impact factors corresponding to 
the overlap between the wave-functions of scattered particle with the wave- 
function of A^— Reggeons, whereas ajy measures the dispersion of the energy 
on the the trajectory around i^h = 0. 

On the other hand, the energy along the trajectories grows with Vh and 
for |fft| — > oo and finally, we have E]\[{ufi;£) ~ Ini/^^. These parts of the 
trajectory give the lowest contribution to the scattering amplitude. 

5.2. Symmetries 

The spectrum of quantized charges q2, ■ ■ ■ ,qN ^s degenerate. This degen- 
eration is caused by two symmetries: 

Qk ^ i-l)^qk (5.4) 

which comes from invariance of the Hamiltonian under mirror permutations 
of particles, H2.3H|1 . and 

qk ^ Qk (5-5) 

which is connected with the symmetry under interchange of the z— and 
z— sectors. Therefore, the four points, {qk}, {{—^)''Qk}, {Qk*} and {(— 1)'^^/;*} 
with k = 2, . . . ,N, are related and all of them satisfy the quantization con- 
ditions (|4.28p and have the same energy. 
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5.3. Descendent states 

Let us first discuss the spectrum along the trajectories with the highest 
conformal charge qn equal zero for arbitrary Uh G M. It turns out [35, 36, 
EZJEj that the wave-functions of these states are built of (A^ — 1)— particle 
states. Moreover, their energies jHH| are also equal to the energy of the 
ancestor {N — 1)— particle states: 

EN{q2,q3,---,qN = 0) = EN.i{q2,q3,---,qN-i)- (5.6) 

Thus, we call them the descendent states of the {N — 1)— particle states. 

Generally, for odd A^, the descendent state '^ ^^~ with the minimal 
energy E]\[{qN = 0) = has for q2 = 0, i.e. for h = 0,1, the remaining 
integrals of motion gs = . . . = ^tv = as well. For /i = 1 + iuh, i.e. q2 / 0, 
the odd conformal charges q2k+i = with k = 1, . . . , {N — l)/2 while the 
even ones q2k / and depend on f^. 

On the other hand, for even A^, the eigenstate with the minimal energy 
^Wiv- ^g |.]^g descendent state of the {N — 1)— particle state which has min- 
imal energy with qN-i + 0. Thus, EJ^''^[qN = 0) = EJ}%{qN-i 7^ 0) > 0. 

Studying more exactly this problem one can obtain ^23 ^| a relation 
between the quasimomenum 9n of the descendent state and the ancestor 
one ^AT-i, which takes the following form 



gifijv 



= -e^^^-i = (-1)^+^ (5.7) 

(jjv=0 



Additionally, one can define a linear operator A jHH E] that maps the 
subspace V^^^ of the (A^—1)— particle ancestor eigenstates with the quasi- 
momentum ^Ar_i = vrA'' into the A''— particle descendent states with qN = 
and 9]\f = 7r(A^ -|- 1) as 

A: y^^_--i-A^) _ vf iv=.(A^+i)) . (5.8) 

It turns out that this operator is nilpotent for the eigenstates which form 
trajectories J2S], i-e. A^^ = 0. Thus, the descendent-state trajectory can 
not be ancestor one for (A^ -|- 1)— particle states. However, it is possible to 
built a single state ^\ with q2 = qs = ■ ■ ■ = qN = 0, i.e. for only one point 
Uh = 0, that has Ej\[ = and the eigenvalue of Baxter Q— operator defined 
as _ 

(5^°(u,u) ~ ^^^, (5.9) 

u 

where a normalization factor was omitted. 

Additional examples of the descendent states will be described later in 

the next sections. 
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6. Quantum numbers of the N = 4 states 

In this Section we present the spectrum for four Reggeons calculated 
by making use of Q— Baxter method. To this end we resum solutions (|4.8p 
and H4.18P numerically, and using them we solve the quantization conditions 
H4.28|l . Earlier, some results for A^ = 4 were only presented in ^\ and some 
numerical results in Ref. ^Hj- Here we show much more data and we present 
more detailed analysis of this spectrum, i.e. resemblant and winding spectra 
of Q3j 94 and corrections to the WKB approximation for A^ = 4. 

For four Reggeons the spectrum of the conformal charges is much more 
complicated than in the three- Reggeon case [20j. Indeed, we have here the 
space of three conformal charges {q2,Q3,Q4)- Thus, apart from the lattice 
structure in q^ we have also respective lattice structures in (73— space. Here 
we consider the case for n^ = so that h = -^ +ii'h- This spectrum includes 
the ground states. For clarity we split these spectra into several parts. We 
perform this separation by considering spectra with different quasimomenta 
9i{q,q) as well as a different quantum number £3, which will be defined in 
the solution H6.2p - (j6.3|l of two quantization conditions H4.28p for A/" = 4 from 
Ref. [IBJ- 

From the first quantization condition of WKB approximation [15j one 
can get the WKB approximation of the charge q^ as 



1/4 _ r^(3/4) 

and the quasimomentum is equal to 



> + 7I^^. 



(6.1) 



9i = -^£ = ^{£2 + 13- h) (mod 27r) , (6.2) 

where i\, ti and £3 are even for even I and odd for odd t. Thus, we have two 
kinds of lattices: with ^4 = 0, vr and with 6*4 = ib7r/2. They are presented in 
Fig. n In these pictures gray lines show the WKB lattice (|6.H1 with vertices 
at £1,^2 £ ^- To find the leading approximation for the charge (73, we apply 
the second relation of WKB approximation jl3| gives 

Im 4^ = (£1-^2-^) =^3- (6.3) 

94 

Notice that the system H6.ip and (|6.3jl is underdetermined and it does not 
fix the charge ^3 completely ^ 



It turns out that after choosing one value of ^4, the lattice in q^ —space 
is still spuriously degenerated^ and this degeneration also corresponds to 

^ degeneration in the leading order of the WKB approximation 
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Figure 1. The spectrum of the integrals of motion 54 for iV = 4 and the total spin 
h — 1/2. The left and right panels correspond to the eigenstates with different 
quasimomenta e*^* = ±1 and ±i, respectively. 

1 /9 ^^^^ 

different lattices in ^3 . The parameter £3 which is defined in (|6.3|1 will be 
used to distinguish these different lattices. 

6.1. Descendent states for N = A 

One can notice that for A'' = 4 and n^ = we have the descendent 
states. They appear in sector with the quasimomentum 64 = vr, which 
agrees with (|5.7|1 . The wave- functions of the descendent states are built of 
three-particle eigenstates with ^3 = 0. Additionally, the spectrum of qs for 
these three-Reggeon states and the descendent state for A/" = 4 is the same 
and it is depicted in Fig. [21 on the left panel. Moreover, the energy of these 
descendent state and the corresponding three-Reggeon states are also the 
same (|5.fi|l . 



6.2. Lattice structure for qs = 

Let us consider the spectrum with 94 7^ and qs = 0, see the right panel 
of Fig 121 In this case the quasimomentum 64 = and the lattice structure 
include vertices that correspond to the ground state. 

Similarly to the A^ = 3 case |2D^ we have in the qj —space a lattice with 
a square-like structure described by (|6.ip . In this case even numbers ii and 
£2 satisfy ^1 -F .£2 G 4Z. Thus, we have the WKB formula 



[qj 



WKB/ 



,1/4 



^Af=4 



2a/2 ^2a/2 



(6.4) 
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Figure 2. The spectra of the conformal charges for iV = 4 and comparison to the 
WKB expansion. On the left panel, the spectrum of q^ with 94 = corresponding 
to the descendent states with 64 = tt. On the right panel, the spectrum of 94 for 
h = 1/2 and qs = with 64 = 0. The WKB lattices are denoted by the gray lines. 



where the vertices are placed outside a disk around the origin of the radius 

-1 



A. 



43/4 



vr 



dx 



Vi 



x" 



r^(3/4) 

2^ 



0.423606 . 



(6.5) 



As before, the leading-order WKB formula (|6.4jl is valid only for \q^ \ ^ 

I 1/2 1 

The energy is lower for points which are nearer to the origin. Similarly 
to the A^ = 3 case, the spectrum is also built of trajectories which extend 
in the (1/^,53, (74)— space. Namely, each point on the qj —lattice belongs to 
one specific trajectory parameterized by the set of integers {£i,i2, ■ ■ ■}■ 

The ground state for A^ = 4 is situated on a trajectory with (^1,^2) = 
(4,0). We find that for this trajectory q^ = Img4 = 0, whereas Re[g'4] and 
E4 vary with u^ as we show in Figure IHl An accumulation of the energy 
levels in the vicinity of i^/j = is described by Eq. (|5.2|l with the dispersion 
parameter CT4 = 5.272. 

On the ^4 —plane the ground state is represented by four points with 
the coordinates (^1,^2) = (i4, 0) and (0,±4). Due to a residual symmetry 
^4 <-^ exp(i/i:7r/2)g4 , they describe a single eigenstate with 



ground 
^3 



0. 



ground 
I4 



0.153589... 



Tjiground 



-2.696640 . . . 



(6.6) 
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0.6- 



0.5- 






0.3 




Figure 3, 
with 92 = 
foriV = 



The dependence of the energy, — -E4/4, and the quantum number, 94/52, 
1/4 + i^,^, on the total spin h = 1/2+ii/h along the ground state trajectory 



with h = 1/2. It has the quasimomentum 64 = and, in contrast to the 
N = 3 case, it is unique. 



(^1/2,^2/2) 


(^ exact) 1/4 


(^WKB)V4 


-^4/4 


(2,0) 


0.626 


0.599 


0.6742 


(2,2) 


0.520 + 0.520i 


0.599 + 0.599 i 


-1.3783 


(3,1) 


0.847 + 0.268 i 


0.899 + 0.299 i 


-1.7919 


(4,0) 


1.158 


1.198 


-2.8356 


(3,3) 


0.860 + 0.860 i 


0.899 + 0.899 i 


-3.1410 


(4,2) 


1.159 + 0.574i 


1.198 + 0.599 i 


-3.3487 



Table 1. Comparison of the exact spectrum of q^ at (73 = and h — 1/2 with the 
approximate WKB expression 1^31 The last column shows the exact energy E4. 

The comparison of H6.4p with the exact results for q^ at h = 1/2 is shown 
in Figure 12 and Tabled One can see that the WKB formula H6.4p describes 
the spectrum with a good accuracy. 



6.3. Resemblant lattices with 







In the previous Section we introduced the parameter £3 which helps us to 
distinguish different lattices. This parameter takes even values for 64 = 0,tt 

and odd ones for 64 = ±7r/2. Let us take £3 = 0. It turns out that in this 

1/2 
case the spectrum lattice in the q^ —space is similar to the corresponding 

lattice in the qj —space, i. e. considering only the leading order of the WKB 
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approximation the q^ —lattice in comparison to the qj —lattice is rescaled 
by some real number. An example of such a lattice for ^4 = is shown 
in Figure 0] One can notice that the non-leading corrections to the WKB 
approximation cause the bending of the lattice structure in Fig. 0] for q^ 

concave whereas for q^ convex. Moreover, we can see the q^ —lattice as 

1/2 
well as ^3 —one have a similar structure to the lattice with 93 = presented 

in Fig. 121 These lattices also do not have the vertices inside the disk at the 

origin Q'a = 54 = 0. 
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Figure 4. The spectra of the conformal charges for A^ = 4 with ^4 = 0, .^3 = 

1 /4 

and ii — 1. On the left panel the spectrum of q^' , while on the right panel the 
spectrum of qj 

Substituting 



Qs = rs e "^3 



and ql^*=rie'^^ 



(6.7) 



into (|6..S|) we obtain a condition for the leading order of the WKB approxi- 
mation 

2 



sin (2(03 - (j)4)). 



(6. 



Thus, for £3 = and for a scale A = r3/r4 > we have (ps = (j)^. This means 

1/2 
that the vertices on the q^ —lattice have the same angular coordinates as 

those from the qj —lattice. Looking at the numerical results in Figure |3] 
we notice that the missing quantization condition for £3 = in the leading 
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Figures. The spectra of the conformal charges for A^ = 4 with ^4 = 0, ^3 = 
and ii = 2. On the left panel the spectrum of qj , while on the right panel the 
spectrum of (73 



1/2 



WKB order should have a form similar to 



93 



1/2 
94 



V4 



(6.9) 



where A^^ = r3/r4 G M is a constant scale for a given lattice ^4. 

It turns out that for a specified quasimomentum we have an infinite 
number of such lattices. They differ from each other by a given scale A^^. 
For example for ^4 = we have another lattice, shown in Fig [3 Its scale 
A2 differs from the scale Ai of the lattice from Fig. 01 The resemblant 
^4 —lattices for ^4 = are described by (jB.ll) with the integer parameters 
ii and £2 satisfying h+h G 4Z. 

The similar lattices also exist in the sector with the quasimomentum 6*4 = 
vr. Some of them are presented in Figures El and For 64 = n the resemblant 
^4 —lattices are also described by H6.H1 but the integer parameters £1 and 
£2 satisfy ^1 -F ^2 S 4Z -F 2. 



6.4- Winding lattices with ^3 7^ 

In the case with £3 7^ we have much more complicated situation than 
for .^3 = 0. According to (|6.8|1 the angles (/>3 and (/>4 defined in H6.7|l are no 
more equal. Moreover, they start to depend on the scale A = r^/r^. 
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Figure 6. The spectra of the conformal charges for TV = 4 with 64 — n, £3 — 
and £4 = 1. On the left panel the spectrum of q/ , while on the right panel the 



spectrum of q^ 
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Figure 7. The spectra of the conformal charges for N — A with 64 — n, £3 — 
and £4 = 2. On the left panel the spectrum of ^4 , while on the right panel the 
spectrum of q^ 



An example of such a lattice, with £3 = 1 and 64 = —n/2 we show in 

FigurelHl For this case, the qj —lattice is defined by l|6.1|l with £1 and £2 odd 
integer numbers satisfying £±+£2 £ 4Z. In Figure |H1 in order to present the 
correspondence between the q^ — and q^ —lattice, we depict only some 
vertices of the lattice which extends in the whole plane of the conformal 
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Figures. The winding spectrum of the conformal charges for TV = 4 and h = 1/2 



1/4 



with 04 = — 7r/2 and £3 = 1. On the left panel the spectrum of q^ , while on the 
right panel the spectrum of q^ 



charges except the place nearby the origin 94 = gs = 0. As we can see the 

1/2 1/2 

gg —lattice is still a square-like one but it winds around the origin gg = 0. 

Looking at Figure [HI let us start from (/)3 = and ^4, = — 7r/4 where 4>^ and 

04 are defined in H6.7|l . Thus, the difference (ps — ip^ = 7r/4 so that our scale 

at the beginning is A = \/2. In this region the vertices of the q^ —lattice are 
in the nearest place to the origin. When we go clockwise around the origin 
of the lattices by decreasing (p^, we notice that ^3 also decreases but much 
slower. Thus, according to H(i.8|l the difference (ps — 1^4 changes. Moreover, 

1 /9 

due to H6.8p . the scale A also continuously rises. Therefore, the q^ —lattice 
winds in a different way than the qj —lattice. After one revolution the 
vertices of the q^ —lattices are at similar places as those with (/>4 decreased 

by 2tt. This provides additional spurious degeneration in q^ . However, 

1/2 
after revolution by the angle 27r, the vertices in q^ —space have completely 

1/2 
different conformal charges qs. The spurious degeneration in the ^3 —lattice 

does not appear. 

We have to add that we obtain the second symmetric structure when 
we go in the opposite direction, i.e. anti-clockwise. Moreover, the winding 
lattice, like all other lattices, extends to infinity on the q^ — and qj —plane 
and does not have vertices in vicinity of the origin, qs = q^ = 0. However, 
the radii of these empty spaces grow with (ps^^. 

Additionally, due to symmetry of the spectrum (|5.5p we have a twin 
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lattice with qk — > ql- Furthermore, the second symmetry (j5.4|l produces 

1/2 

another lattice rotated in the ^3 —space by an angle 7r/2. Notice that this 
symmetry H5.4p exchanges quasimomentum ^4 <-> 27r — 04(mod 27r). Thus, 

the spectrum for 9^ = 7r/2 is congruent with the spectrum of 6*4 = — 7r/2 but 

1/2 
it is rotated in the q^ —space by 7r/2. 
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Figure 9. The winding spectrum of the conformal charges for TV = 4 with h — 1/2, 
64 = and £3 = 2. On the left panel the spectrum of q^' while on the right panel 
the spectrum of ^3 

As we said before for 9^ ± 7r/2 we have winding spectra with odd £3. 
Similarly, for even ^3 7^ we have winding spectra with 9^ = 0, vr, thus, 
the lattice with the lowest non-zero is corresponds to l^sl =2. We present 
some points of this spectrum in Fig. El for which ii, £2 in (|6T|) are even and 
^1 + ^2 £ 4Z + 2. In this case we start with 03 = 7r/4 and (p^ = which 
implies the beginning scale A = \/2. The spectra wind as in the previous 
case. 

Similarly, for 64 = vr we have also spectrum with the lowest l^sl = 2. 
It is defined by (|f).l|) with ii, £2 even and ^1+^2 S 4Z. In this case the 
angles (/>3 = and 4)4 = —tt/A so we also have the beginning scale A = \/2, 
as depicted in Figure ^1 In order to describe the winding spectra better 
we may introduce an integer parameter £4 which helps us to number the 
overlapping winding planes of the spectra and name spuriously-degenerated 
vertices in the q^ —plane. 

To sum up, even for a given quasimomentum we have many lattices which 

1/2 
overlap, so that vertices of the lattices, especially in ^3 —space, make an 

impression of being randomly distributed. However, as we have shown above. 
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Figure 10. The winding spectrum of the conformal charges for iV = 4 with h — 1/2, 
04 = TT and l^ = 2. On the left panel the spectrum of qj while on the right panel 



the spectrum of q^ 



1/2 



those spectra may be distinguished and finally described by (|6.H1 and H6.3J1 . 
Still, there is a lack of one nontrivial WKB condition which would uniquely 
explain the structure of the resemblant and winding lattices. 



6.5. Corrections to WKB 
Let us consider the spectrum of the conformal charge q^ iov N = A with 



93 



and h 



l+Uh 



2 — . It turns out that for n/j 7^ it has similar square- 
like lattice structure like that with n^ = 0, see Fig. El on the right panel. 
Similarly to the case with three reggeized gluons [201 ^^ have evaluated the 
conformal charges 54 with even n/j with high precision. We have done it 
separately for q^ with Im[g4] = and Re[g4] = 0. Next, we have fitted 



(r) 



,« 



coefficient expansion of the WKB series, a), and a\ , respectively. 
In ^ni the series formula for q^ looks as follows 



1/4 
94 



TT' 



3/2 



2r2(i/4) 



Q{n) 



1 + 



\Q{n) 



+ Y,^k 



k=2 



IQin) 



where 



j7T{2k~l)/4 

and £1,^2, n = {ni, . . . ,nj^} are integer. 



fc=i 



— + i— 



(6.10) 



(6.11) 
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nh 


coef. 


k = 2 


k = 3 


k = 4 


k = 5 







2.9910566246 


-24.021689 


91.591 


645.5 




4- 


-4.9910566246 


28.021689 


-148.830 


1656.7 




«r+4- 


-2.0000000000 


4.000000 


-57.239 


2302.2 


2 


4" 

4" 


-1.3991056625 


4.674008 


-4.516 


-95.7 




-0.6008943375 


-0.674008 


20.388 


-200.8 




4" +4" 


-2.0000000000 


4.000000 


15.872 


-296.5 


4 


4" 


-1.351212983 


3.462323 


-11.322 


43.164 




-0.648787017 


0.537677 


-0.096 


0.611 




4''+4' 


-2.000000000 


4.000000 


-11.418 


43.775 


6 


4" 


-0.8882504145 


1.883461 


-5.4248 


18.081 




4' 


-1.1117495855 


2.116539 


-4.8117 


12.091 




4" +4" 


-2.0000000000 


4.000000 


-10.2365 


30.172 


8 


4" 


-0.6326570719 


1.197694 


-2.94579 


8.3372 




-1.3673429281 


2.802305 


-5.96129 


12.1238 




4''+4' 


-2.0000000000 


4.000000 


-8.90708 


20.4610 



Table 2. The fitted coefficient to the series formula of q^ Il6.1fl|l with Uh — 0, 2, 4, 6 
and 8 



Here for N 
different expansion parameter 



4, similarly to the A^ = 3 case in Ref. J2D], we have a 
|Q(n)| and in this case the parameter 



(6.12) 



is decreased by 5/4. The expansion coefficients of the series (j6.1()p are shown 
in Table |21 Similarly to ^H] the coefficients ao = 1 and ai = 1. The 
remaining coefficients depend on n^. Moreover, the coefficients o^ with 
A: > 1 are different for real q^ and for imaginary qj but one may notice 



At) 



,W 



that for A; = 1,2 the sum al. + a^. G Z and does not depend on rih. Thus, 
to describe the quantized values of q^ more generally one has to use both 

(r) (i) 

sets of the coefficients, a^ and o^ , or perform the expansion with two small 
independent parameters, i.e. g2/l2('^)P ^^'^ l/IQl*^)!^- 

Using the series (|6.1()|1 with (|6.12|1 and coefficients from Table [21 gives 
good approximation of the conformal charges q^ with q^ = 0. However, in 
order to have a better precision one has to introduce an additional expansion 
parameter. 
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7. Summary and Conclusions 

In this paper we discussed four Reggeon states which appear is scattering 
amphtude of strongly interacting particles in high energy Regge limit of 
QCD (|2.ip . i.e. in Generalized Leading Logarithm Approximation (GLLA) 
[HlEBISni- Due to the colour factor description of A^— Reggeon system is 
a complicated task especially for A^ > 3. To this end multi-colour limit 
is performed. In this limit the Reggeon Hamiltonian separates into two 
one-dimensional equations and as a result the multi-Reggeon system get 
completely solvable. However, in order to solve it one has to use advanced 
integrable-model methods. 

In the multi-colour limit the equation for the A^— Reggeon wave-function 
takes a form of Schrodinger equation H2.1()p for the non-compact XXX Heisen- 
berg magnet model of SL(2, C) spins s ^OlESl^B- Its Hamiltonian describes 
the nearest neighbour interaction of the Reggeons 02 ESI propagating in the 
two-dimensional transverse-coordinates space (|2.6|1 . The system has a hid- 
den cyclic and mirror permutation symmetry H2.36J) . It also possesses the 
set of the (A^ — 1) of integrals of motion, which are eigenvalues of conformal 
charges ^2j, q^ and q^, H2.H3|1 - I|2.H4|1 . Therefore, the operators of confor- 
mal charges commute with each other and with the Hamiltonian and they 
possess a common set of the eigenstates. 

To solve the Reggeon problem for more than three particles one uses the 
more sophisticated technique, i.e. the Baxter Q— operator method j2Zj. It re- 
lies on the existence of the operator Q{u,u) depending on the pair of complex 
spectral parameters u and u. The Baxter Q— operator has to commute with 
itself (|3.H1 and with the conformal charges (|H.2jl . It also has to satisfy the 
Baxter equations H,S.,Sp - (|,S.4|l . Furthermore, the Q— operator has well defined 



analytical properties, i.e. known pole structure H3.12|l and asymptotic be- 
haviour at infinity H3.14|l . The above conditions fix the Q— operator uniquely 
and allow to quantize the integrals qk- In turns out J22] that the Reggeon 
Hamiltonian can be rewritten in terms of Baxter Q— operator (j3.15p . There- 
fore, combining together the solutions of the Baxter equations and the con- 
ditions for qk with the Schrodinger equation, we can calculate the energy 
spectrum H3.24p . Moreover, we are able to determine the quasimomentum 
of the eigenstates (|3.2()p . i.e. the observable which defines the properties of 
the state with respect to the cyclic permutation symmetry (|2.36p . 

In order to obtain the exact values of qk we have used the method [T5|[T7] 
which consists in rewriting the Baxter equations and the other conditions for 
Q— operator eigenvalues as the A^— order differential equation (|4.2p . Solving 
this equation one obtains the conditions for quantized qk l|4.28p and the 
formulae for the energy H4.15p and the quasimomentum H4.25p . 

In the case with A^ = 4 particles we have constructed the spectrum for 
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Uh = depicting complicated interplay between the conformal charges: q^ 
and Qi. Such a complete analysis has been performed here for the first time. 
Earlier, for A^ = 4 full spectrum of q^ was shown in Ref. ^21 (however the 
corresponding gs spectrum was not discussed) and some values of q^ were 

found in Ref. jTHI. The spectrum of ^4 has a structure of square-like 
lattice (|().4|1 . Fig. ^ In this case the spectrum of qs is more complicated. 
It turns out that it has a few possible forms. Firstly, there are simple 
states with q^ = 0, Fig. |2j Moreover, we have found the q^ —lattices 
whose distribution of vertices is similar to the distribution of vertices of 
^4 —lattice. We have called this structures as resemblant lattices. Figs. ^ 

1 /9 

13 Secondly we have demonstrated the existence of q-^ —lattices that may 
be called winding lattices. They wind around the origin and in course of 
this winding the distance between vertices increases and also the lattice goes 
away from the origin, Figs. lSllTIll The WKB approximation does not describe 
these structures (|H.:-{|) exactly because one quantization condition is missing. 

Finally, we have also considered the descendent states with q^ = 0, Fig. |21 

1/3 
They have quasimomentum 84 = it and the structure of their ^3 —lattices is 

1 /s 

the same as for q^ —lattices in the three- Reggeon case with 64 = [^1^ . 

For the A^ = 4 Reggeon states with q^ = we have calculated corrections 
to the WKB approximation H6.1()|l , Table El Similarly to the A'^ = 3 case 
|2()| they come from the fact of using only one expansion parameter r/ for 
two different conformal charges, q^ and ^2- One has to notice that these 
corrections are subleading in comparison to the WKB limit |15| . 

The Q— Baxter method is complicated. However, it allows us to solve 
the reggeized gluon state problem for more than N = S particles. The 
above calculations are of interest not only for perturbative QCD but also 
to statistical physics as the SL(2,C) non-compact XXX Heisenberg spin 
magnet model (13 ESI Ell- 
in order to find the full scattering amplitude one has to calculate the 
impact factors, i.e. overlaps between the Reggeon wave function and the 
wave functions of the scattered particles. These impact factors strongly 
depend on the scattered system and they are hard to calculate due to a large 
number of integrations. Thus, this work can be a first step in computing of 
contribution of the four-Reggeon state to the hadron scattering amplitude. 

Acknowledgements 

I would like to warmly thank to Michal Praszalowicz for fruitful discus- 
sions and help during writing this work. I am very grateful to G.P. Korchem- 
sky, A.N.Manashov and S.E. Derkachov whom I worked in an early state of 
this project. I also thank to Jacek Wosiek for illuminating discussions. This 



phdn4a printed on 1st February 2008 39 



work was supported by KBN PB 2-P03B-43-24, KBN PB 0349-P03-2004-27 
and KBN PB P03B-024-27:2004-2007. 



References 

[1] M. Gell-Mann, M.L. Goldberger, F.E. Low, E. Marx, F.Zachariasen. Elemen- 
tary particles of Conventional Field Theory as Regge poles. III. Phys. Rev. B, 
133:145, 1964. 

[2] M. Gell-Mann, M.L. Goldberger, F.E. Low, E. Marx, F.Zachariasen. Elemen- 
tary particles of Conventional Field Theory as Regge poles. IV. Phys. Rev. B, 
133:161, 1964. 

[3] V. N. Gribov. A Reggeon diagram technique. Sov. Phys. JETP, 26:414-422, 
1968. 

[4] V. S. Fadin, E. A. Kuraev, and L. N. Lipatov. On the Pomeranchuk singularity 
in asymptotically free theories. Phys. Lett., B60:50-52, 1975. 

[5] J. Bartels. High-Energy behavior in a nonabelian gauge field theory. Phys. 
Lett, B68:258, 1977. 

[6] H. Cheng and C. Y. Lo. High-energy amplitudes of Yang-Mills theory in 
arbitrary perturbative orders. 1. Phys. Rev., D15:2959, 1977. 

[7] L. N. Lipatov I. I. Balitsky. The Pomeranchuk singularity in Quantum Chro- 
modynamics. Sov. J. Nucl. Phys., 28:822-829, 1978. 

[8] E. A. Kuraev, L. N. Lipatov, and V. S. Fadin. The Pomeranchuk singularity 
in nonabehan gauge theories. Sov. Phys. JETP, 45, 1977. 

[9] J. Bartels. High-energy behavior in a nonabelian guage theory. 2. First cor- 
rections to r(ri — > to) beyond the leading Ins approximation. Nucl. Phys., 
B175:365, 1980. 

[10] J. Kwiecinski and M. Praszalowicz. Three gluon integral equation and odd C 
singlet Regge singularities in QCD. Phys. Lett., B94:413, 1980. 

[11] T. Jaroszewicz. High-energy multi - gluon exchange amplitudes. Triest 
preprint IC/80/175. 

[12] R. A. Janik and J. Wosiek. Solution of the odderon problem. Phys. Rev. Lett., 
82:1092-1095, 1999. 

[13] P. Gauron, B. Nicolescu, and L. Szymanowski. A possible field theoretical 
description of the odderon. IPNO/TH 87-53. 

[14] L. Lukaszuk and B. Nicolescu. A possible interpretation of p p rising total 
cross- sections. Nuovo Cim. Lett., 8:405-413, 1973. 

[15] S. E. Derkachov, G. P. Korchemsky, and A. N. Manashov. Noncompact Heisen- 
berg spin magnets from high-energy QCD. Ill: Quasiclassical approach. Nucl. 
Phys., B661:533-576, 2003. 

[16] G. P. Korchemsky, J. Kotanski, and A. N. Manashov. Compound states of 
reggeized gluons in multi-colour QCD as ground states of noncompact Heisen- 
berg magnet. Phys. Rev. Lett., 88:122002, 2002. 



40 phdn4a printed on 1st February 2008 



[17] S. E. Derkachov, G. P. Korchemsky, J. Kotanski, and A. N. Manashov. Non- 
compact Heisenberg spin magnets from high-energy QCD. II: Quantization 
conditions and energy spectrum. Nud. Phys., B645:237-297, 2002. 

[18] H. J. De Vega and L. N. Lipatov. Exact resolution of the Baxter equation for 
reggeized gluon interactions. Phys. Rev., D66:074013, 2002. 

[19] H. J. De Vega and L. N. Lipatov. Interaction of reggeized gluons in the 
Baxter-Sklyanin representation. Phys. Rev., D64:114019, 2001. 

[20] J. K otanski. Three particle Pomeron and odderon states in QCD. 2005. 
piep-th/0603238, 

[21] Jan Kotanski. Reggeized gluon states in Quantum Chromodynamics. 2005. 
The PhD Thesis, hep-th/0511279 

[22] M. A. Braun. The interaction of reggeized gluons and Lipatov 's hard pomeron. 
- lectures. University of Santiago de Compostela, 17506 Santiago de Com- 
postela, Spain, 1994. 

[23] G. 't Hooft. A planar diagram theory for strong interactions. Nucl. Phys., 
B72:461, 1974. 

[24] L. N. Lipatov. Pomeron and odderon in QCD and a two-dimensional conformal 
field theory. Phys. Lett., B251:284-287, 1990. 

[25] L. N. Lipatov. High-energy asymptotics of multicolor QCD and two- dimen- 
sional conformal field theories. Phys. Lett., B309:394-396, 1993. 

[26] R. J. Baxter. Exactly Solved Models in Statistical Mechanics. Academimc 
Press, London, 1982. 

[27] S. E. Derkachov, G. P. Korchemsky, and A. N. Manashov. Noncompact Heisen- 
berg spin magnets from high-energy QCD. I: Baxter Q-operator and separation 
of variables. Nucl. Phys., B617:375-440, 2001. 

[28] P. Di Francesco, P. Mathieu, D. Senechal. Conformal Field Theory. Springer, 
New York, 1982. 

[29] J. B. Zuber. An introduction to Conformal Field Theory. Acta Phys. Polon., 
B26:1785-1813, 1995. 

[30] E. K. Sklyanin. Quantum inverse scattering method. Selected topics. 1991. 



hep-th/9211111 



[31] L. D. Faddeev. Algebraic aspects of Bethe Ansatz. Int. J. Mod. Phys., 
A10:1845-1878, 1995. 

[32] L. D. Faddeev. How Algebraic Bethe Ansatz works for integrable model. 1996. 



hep-th/9605187 



[33] L. D. Faddeev, E. K. Sklyanin, and L. A. Takhtajan. The quantum inverse 
problem method. 1. Theor. Math. Phys., 40:688-706, 1980. 

[34] J. Kotanski and M. Praszalowicz. Solutions of the quantization conditions for 
the odderon charge q^ and conformal weight h. Acta Phys. Polon., B33:657- 
682, 2002. 

[35] G. P. Vacca. Properties of a family of n reggeized gluon states in multicolour 
QCD. Phys. Lett., B489:337-344, 2000. 



phdn4a printed on 1st February 2008 41 



[36] L. N. Lipatov. Duality symmetry of Reggeon interactions in multicolour QCD. 
Nucl Phys., B548:328-362, 1999. 

[37] J. Bartels, L. N. Lipatov, and G. P. Vacca. A new odderon solution in pertur- 
bative QCD. Phys. Lett, B477:178-186, 2000. 

[38] G. P. Korchemsky. Bethe ansatz for QCD pomeron. Nucl. Phys., B443:255- 
304, 1995. 

[39] T. Jaroszewicz. Infrared divergences and Regge behavior in QCD. Acta Phys. 
Polon., Bll:965, 1980. 

[40] L. A. Takhtajan and L. D. Faddeev. The quantum method of the inverse 
problem and the heisenberg xyz model. Russ. Math. Surveys, 34:11-68, 1979. 

[41] N. M. Bogolyubov, A. G. Izergin, and V. E. Korepin. Quantum inverse scat- 
tering method and correlation functions. Univ. Press, Cambridge, 1993. 

[42] L. N. Lipatov. High-energy asymptotics of multicolor qcd and exactly solvable 
lattice models. JETP Lett., 59:596-599, 1994. 

[43] L. D. Faddeev and G. P. Korchemsky. High-energy QCD as a completely 
integrable model. Phys. Lett., B342:311-322, 1995. 



